Geometry of lipid vesicle adhesion 
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The adhesion of a lipid membrane vesicle to a fixed substrate is examined from a geometrical point 
of view. This vesicle is described by the Helfrich hamiltonian quadratic in the mean curvature; it 
interacts by contact with the substrate, with an interaction energy proportional to the area of 
contact. We identify the constraints on the geometry at the boundary of the shared surface. The 
result is interpreted in terms of the balance of the force normal to this boundary. No assumptions 
are made either on the symmetry of the vesicle or on that of the substrate. The strong bonding 
limit as well as the effect of curvature asymmetry on the boundary are discussed. 

PACS numbers: 87.16.Dg, 46.70.Hg 
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I. INTRODUCTION 

Geometrical models provide a surprisingly robust phe- 
nomenological description of the equilibria of physical 
membranes 0, |[ ^ ^ . The hamiltonian which de- 
scribes the membrane is constructed as a sum of geomet- 
rical scalars; in particular, lipid membrane vesicles are 
well approximated by one which is quadratic in the mean 
curvature §, 0, |. Such models can also be extended to 
model adhesion between vesicles or between a vesicle and 
a rigid substrate ^ |ll|, processes which are increas- 
ingly relevant to biophysics. (Two reviews are jlj, ^.) 
In this paper, we examine one important aspect of this 
problem, the geometry of the contact boundary, which, 
surprisingly, does not appear to have been examined in 
any generality. 

To model the interaction one can exploit, as for an iso- 
lated membrane, the geometrical scalars characterizing 
the surface of contact as well, perhaps, as its boundary. 
As such, this is not a model of the adhesion of individual 
molecules to specific sites on the membrane, a task which 
lies beyond the scope of this continuum description. 

In its simplest form, which is the one we consider, the 
interaction hamiltonian is proportional to (minus) the 
area of contact. The energy associated with the bound- 
ary of the contact region is ignored. Axially symmet- 
ric configurations have been studied thoroughly in this 
'ideal' context||ll|, Q . In the adhesion of 'linear' 
vesicles in two dimensions was considered. More recently, 
and 
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p8| perturbation theory has been devel- 
oped in the strong bonding limit, in which the bending 
energy itself is small compared to that of adhesion. Non- 
axially symmetric configurations of an adhering vesicle 
under the effect of gravity were studied in |ll9| using nu- 
merical techniques. We note that a more realistic treat- 



ment of adhesion considering chemically structured or 
rough surfaces has been provided in pc[ |. 

For definiteness, we will assume that one of the inter- 
acting surfaces is a fixed substrate, although it is simple 
to relax this restriction. We will not, however, assume 
that this surface is flat. More significantly, we relax the 
assumption that the vesicle geometry is axially symmet- 
ric. In previous axially symmetric work the geometric 
origins of the boundary conditions are not clear, because 
the technique is tailored so finely to the symmetry; nor 
is it clear to what extent they will survive the relaxation 
of symmetry. This is a less than desirable situation in a 
model which is intrinsically geometric to begin with. Of 
course, one is also interested in geometries which are not 
axially symmetric: to mention just one context where 
this would be the case, we note that all configurations 
with a negative area difference appear to be inconsistent 
with axial symmetry plf . Indeed, it may also be energet- 
ically favorable for an initially axially symmetric vesicle 
to adhere to a substrate in a manner which breaks its 
original symmetry. 

Our first approach will be to search for minima of 
the energy. To do this, we will exploit the geomet- 
rical framework introduced recently to describe lipid 
membranes [p^, and extended to accomodate edge effects 
in . The extension to adhering geometries introduces 
its own subtleties due to potential discontinuities at the 
boundary of contact: the energy is stationary only when 
appropriate constraints on the vesicle geometry are sat- 
isfied on this boundary. Our treatment of the problem is 
divided into three parts. To establish our bearings, we 
begm m Sect. |l| with a rederivation of the Young equa- 
tion for a liquid droplet where the bonding to the wall 
competes with the surface tension of the drop. In Sect. 
Ill, we consider lipid vesicles described by the Helfrich 
hamiltonian. Discontinuities at the boundary of the con- 
tact region are discussed in Sect. [II A. In Sect. IIIB, 
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tablishing contact with the recent work in this direction, 
we study the strong bonding limit in which the bend- 
ing energy is ignored but any asymmetry between the 
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layers is accounted for. In this limit, discontinuities at 
the boundary of the contact region are not smoothed; a 
non-vanishing contact angle does not imply a divergent 
energy. Finally, in Sect. IV we consider the general case. 
The finiteness of the curvature energy necessarily elim- 
inates an angle of contact between the vesicle and the 
substrate; stationary energy completely fixes the curva- 
ture at the boundary. This generalizes the situation for 
axially symmetric shapes, where, as is well known, the 
curvature is completely fixed at the boundary the 
vesicle radius of curvature normal to the line of contact 
is inversely proportional to the square root of the bond 
strength, the tangential radius as well as the potentially 
non- vanishing off-diagonal curvature are both completely 
fixed by the substrate geometry. We demonstrate that 
the boundary condition is not modified by a curvature 
asymmetry. In Sect. we use the expressions for the 
internal stress tensor in a lipid membrane developed in 
p2| , to provide a surprisingly simple interpretation of the 
boundary condition in terms of the balance of forces at 
the boundary for a flat substrate. We end with some 
brief conclusions. 



II. SURFACE TENSION DOMINATED MODEL 

It is worthwhile to first review the adhesion of a drop of 
liquid of fixed volume onto a surface. Here the focus is on 
the competition between the surface tension of the liquid 
and the attraction between the liquid surface and the 
substrate. The former tends to reduce the surface area 
of the drop; the latter to increase the area of contact. 
The energy is given as a sum of three terms 



F ^ ^J.A~ W^contact " piV ~ Vq) . 



(1) 



The energy associated with the constant surface tension 
/i of the drop is proportional to its total surface area A; 
that associated with adhesion is proportional to (minus) 
the area of contact ^contact between the drop and the 
substrate. The parameter w is the attractive contact po- 
tential. The third term involving the Lagrange multiplier 
p implements the volume constraint fixing the enclosed 
drop volume V at the value Vq. 

The equilibrium drop configurations are those at which 
the energy (|I|) is stationary. The problem posed here dif- 
fers from the standard isoperimetric problem in that the 
area at different positions on the surface gets weighted 
according to whether or not it lies in the contact region, 
which itself is determined by the outcome of the vari- 
ational problem. Indeed the contact surface might be 
weighted negatively. For physically realistic parameters, 
however, an equilibrium is realized. The energy is always 
bounded from below. 

In equilibrium, the curvature of the drop's surface will 
suffer a discontinuity along the boundary C of the con- 
tact region. We parameterize the embedding of the free 
surface of the droplet in three dimensional space as fol- 
lows: X = X(^°), and the substrate x = X(^"), a = 1, 2. 




x=X( ^ , 4 ) 



FIG. 1: Definition of the quantities used in the description of 
the geometry of adhesion. 



The energy is a functional both of X for the free surface 
and X for the region of contact. They coincide on C, 
X = X. See Fig. 1. We now recast the first two terms 
appearing in F as (Aftoc is the area of the free surface) 



ly^contact = A*^frco + (M " t«)^contact 

-'2'' ^ ■ ' ^ f -^2^ ^2) 



free J contact 
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Here 7 is the determinant of the metric ^ab induced on 
the free surface given by jab = ■ ^b, where = cJ^X 
are vectors tangent to the surface Similar definitions hold 
for the geometrical quantities, indicated with a tilde, as- 
sociated with the substrate. Note that the boundary C 
may possess disconnected components. 

To derive the equations describing the equilibrium 
shape of a droplet, let us consider a variation of the 
embedding of the free surface, X ^ X + (5X. We let 
n denote the unit normal to the free surface. We de- 
compose the displacement with respect to the spatial ba- 
sis adapted to this surface, {ea,n}: 5X = $°ea -I- <i>n. 
We have for the corresponding variation of the induced 
metric Sx-fab = Sifab^" + Va^b + Vb$a- The normal de- 
formation is proportional to the extrinsic curvature ten- 
sor, Kab = Gb • dan. The mean extrinsic curvature is 
K = Kabj""^- The tangential deformation is the Lie 
derivative of jab along the vector field $"; Vq is the co- 
variant derivative compatible with jab- 

On the boundary C, the fixed substrate constrains i5X 
to lie along the contact region. We will ignore this for the 
moment, treating (SX as though it were unconstrained on 
C. Then variation of ^ftoc gives. 



SxAfj^cc — 



dAK<^ 



free 



ds r $a . 



(3) 



Here l" is the outward pointing normal to C on the free 
surface; s is arclength along C. We also have that the 
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variation of the enclosed volume is 



6xV 



(4) 



free 



Remote from C only the normal projection of the vari- 
ation plays a role in determining the equilibria of the 
droplet. This is generally true regardless of the model. 

For this model, the free surface satisfies = p, as 
follows from the first term in Eq. (^), together with Eq. 

. Note that there is no boundary term associated with 
the normal deformation (f>. This contrasts with the tan- 
gential deformation whose only net physical effect is to 
induce a displacement of the boundary. 

The boundary deformation we have described is not 
free: the variation JX on C is constrained to lie tangent 
to the substrate. Without loss of generality, we can also 
assume that it is normal to the boundary C, so that 



(5X = $0 z , 



(5) 



where z is the outward unit normal to C on the substrate 
(see Fig. 1). We then have for the integrand appearing 
in the boundary term in Eq.(^), Z^^^ = 1 • 5X = 1 • z $o, 
where 1 = P Gq is the surface vector 1°- treated as a spatial 
vector. The boundary contribution to the variation of the 
free surface ^x^free is then 



ds 1 • z <1> 



• 



(6) 



We now consider the variation of the area of contact 
^eontact- The dcformatiou 5X of the free surface will 
induce a variation in ^contact through the boundary that 
they share, 



(^X^contaet 



(7) 



which is a two dimensional analog of Eq. (Q) . Note that 
the substrate need not be planar. We can now read off 
the total boundary contribution to the variation SxF, 
with F given by Eqs. (|l|), (H). In equilibrium, we require 
that 



/ ds [^1 • z + (^ - $0 = 
Jc 



(8) 



for an arbitrary <i>o. We therefore conclude that 

w = ^(l + l-z). (9) 

Defining the contact angle 8 by cos 8 = — 1 • z, this ex- 
pression reproduces the well known Young equation. 

III. LIPID VESICLE ADHESION 

A lipid membrane is modeled by the Helfrich bending 
energy 



(10) 



For definiteness, we will focus on either of two variants of 
the model: in both versions the enclosed volume and the 
total surface area are fixed; in the spontaneous curvature 
model, K is replaced hy K — Kq in Eq.(^o|) where the 
constant Kq is the spontaneous curvature; in the bilayer 
couple model, the area difference, proportional to the 
integrated mean curvature 



M ^ dAK , 



(11) 



is also fixed ||2l|l. Thus we construct the constrained en- 
ergy, 

F - Fb~wA,onte.ct + ^^{A-Ao)+(3{M-Mo)-p{V-Vo) . 

(12) 

Here fi is now the Lagrange multiplier associated with 
the area constraint; likewise /3 is that associated with 
the area difference constraint. 



A. Discontinuities resolved 

In the simple model discussed in Sect. || there is no 
energy penalty associated with discontinuities along the 
contact boundary C. However, both the intrinsic and 
extrinsic curvature will suffer a discontinuity along C. 
When the curvature of the vesicle contributes to its en- 
ergy, such a discontinuity will generally result in a singu- 
larity in the energy. Because this singular contribution 
has support on C, it is no longer valid to decompose the 
energy into two parts, F = Fftoe + -Feontact, 

This point is well illustrated by considering an axi- 
ally symmetric surface. Cylindrical polar cooordinates 
{p, z,ip} are introduced on R^; constant (p curves on the 
surface are parametrized by arclength i. The surface is 
then described completely once p = R{£) is specified. The 
extrinsic curvature tensor consistent with axial symmetry 



(13) 



where the scalars Ki and Kji are the two principal cur- 
vatures, and 1°" is the outward pointing unit normal to 
the circle of fixed £, £"' — (1,0). We identify the scalar 
curvature TZ = 2detK = 2KiKR, and K = Ki + Kr. 
Now let 9 be the angle which the tangent to a curve of 
fixed if makes with the positive x axis. The principal 
curvatures are then given by 



Kf = 0' 



sm9 



(14) 



The prime denotes a derivative with respect to i. We 
have for the integrated mean curvature 



M = 2tt / dm 9' 



sin 9 



(15) 



Suppose that 9 suffers a discontinuity 8 on the circle at 
i = £o, so that 9{£) « eH{£ - 4), where H is the step 
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function. There is a finite contribution from this circle 
given by 

/ dAK = 2n de9'R = 271 R{eo)Q- (16) 

Jc Jto-f- 

The mean curvature is thus integrable across the discon- 
tinuity. In general we have the decomposition 



M ■ 



I dA K ~ ds arccos(l • z) , 

^drop/C Jc 



(17) 



with cos O ~ —I ■ z, and where the notation for the nor- 
mals is that introduced in Sectjl|. 

We note that the Gauss-Bonnet invariant for a vesicle 
of spherical topology can likewise be decomposed 



drop/C 



dAn + 2 / dsAK = 8TT, 



(18) 



c 



where Ak = k— k is the discontinuity in the geodesic cur- 
vature of C. For an axially symmetric vesicle, the value of 
the latter term is 2tt cos Q. (An analogous decomposition 
of the Hilbert-Einstein action arises in the study of the 
dynamics of thin-shells in general relativity (24j.) As long 
as the adhering vesicle remains intact, the Gauss-Bonnet 
invariant will not play a role in adhesion. Though each 
of the two components appearing in Eq.([l^) will behave 
non-trivially under deformation of the surface, their sum 
will not change. 

The geometric invariant does possess a singular- 
ity at a curvature discontinuity. We note that this sin- 
gularity is identical to that arising from the alternative 
quadratic invariant, KabK"-^. This is because the Gauss- 
Codazzi equation, TZ — — KabK°'^ , identifies their 
difference as the scalar curvature TZ which according to 
Eq.(|l8|) picks up a finite contribution at a discontinuity. 
In an axially symmetric geometry the troublesome term 
in the bending energy is 9'^: 



dAK^ « 2n / diRe'^ + 



(19) 



The 0'^ term appearing in the integrand gives rise to a 
delta function squared singularity across the boundary. 
To eliminate the corresponding divergence in the energy, 
we do require = 0. The surface must be differentiable 
across C. It is straightforward to bootstrap this axially 
symmetric analysis to the general case by introducing 
Gaussian normal coordinates adapted to the boundary. 
In general, we require that 8 = or z = —1. 



B. Strong Bonding Limit 

Before adressing the full problem, let us consider the 
strong bonding limit, a << wA. At lowest order the 
bending energy Ft, is ignored in Eq. (^2|), and the vari- 
ational problem reduces to the the minimization of the 



contact energy subject to the three constraints. (In this 
section, we will use the language appropriate to the bi- 
layer couple model.) Whereas Q necessarily vanishes 
on the contact boundary for the Helfrich hamiltonian, 
it need not vanish in this limit. 

Let us first consider the variational problem on the free 
surface of the vesicle. Under a tangential deformation of 
this surface any scalar function J-', and in particular J- = 
fi + f3K, transforms as a divergence which is transferred 
to the boundary: 



dAT 



[ dsla^'^J'. 

Jc 



(20) 



This is because (5||JF = ^"^da^- The details of are 
irrelevant. Since is constrained to lie tangent to the 
contact region, from Eq. (|^), with J- = (3 we have 
then 



^11 



dA in + PK) 



free 



ds{^i + l3K)\-z'S>o. (21) 



Note that this expression reduces to Eq.(g) when /i = 
1,/? = 0. 

We now examine a normal deformation of the free sur- 
face. The Euler-Lagrange equation which determines the 
local vesicle equilibrium shape is obtained by demanding 
that the energy be stationary with respect to normal de- 
formations of the free surface. Its derivation within the 
pre sent framework has been discussed in detail elsewhere 
|22| . Let us focus on the normal deformation of the new 
ingredient with respect to the model discussed in Sect. || 
appearing in Eq.(|lJ) which is M. We consider the con- 
tributions from the free surface, C and the contact region 
as given by Eq.(p^) separately. We begin with Mfree- We 
have 



S I M, 



free 



dAn<i> 



ds V I $ . 



(22) 



We have used Sj_K = -V^$ - K"-^Kab^, as well as the 
Gauss-Codazzi equation, and — 1°"^ a denotes the 
derivative along 1. It is now straightforward to read off 
the bulk Euler-Lagrange equation. 



(23) 



Note that this equation is second order in derivatives. 

To proceed with the determination of the boundary 
conditions, we need to identify the independent uncon- 
strained variation at the interface. We identify these as 
^0 — z ■ (5X and its derivative along z, Vj_<i>o. We will, 
however, continue to use V_l<I>o to denote 1 • zV_l<I>o. We 
note that the normal deformation at the boundary, using 
Eq. (I) is 

$ = n-(5X = n-z$o, (24) 
and on the boundary C, its normal derivative is 

V_L$ V_L (n • z$o) 

= (n • z) V_L<i>o + ^0 z ■ V_Ln + $o n • V^z 
= (n • i)Vi$o + (1 • i)[Ki_ + (1 • z)ifi]$o(25) 
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where we have defined K± = Katl'^l , and K± = 
KabZ°'z''. We have used the fact that n- Vj^z — {l-z)'^K± 
as well as Vj^n = I'^Ka^eb. Therefore the boundary con- 
tribution to Eq. ([2^ ) takes the form 

5±Affrcc = - / rfs{(n-i)Vi$o+(l-z)[if± + (l-z)i^±]$o}. 
Jc 

(26) 

For the boundary contribution Mc, we have 



IV. 



NO APPROXIMATIONS 



SxMc = / ds [e + 5xQ] 
Jc 



(27) 



We emphasize that this term contributes only to the 
strong bonding limit. The first term comes from the vari- 
ation of arclcngth: 5xds = k^q. We now exploit the fact 
that cos & ~ — 1 • z to express Sx 6 in terms of (5x (1 • z) ; to 
evaluate the latter we note that SxQ = {smQ)^^SxO--z), 
and that 

Sxil-i) = {Sxl)-i^n-iin-Sxl). (28) 

In general, we have for the tangent vectors to the free 
surface, 



n • SxBa = Va^ - i^ab*'' ■ 



(29) 



so that 



(<5jfl)-z= (Vi$-(l-z)X^$o)n-z. (30) 



Note that SxMc need not vanish even when Mc itself 
does. 

Finally, the induced change in M on the contact region 
due to the displacement of C is just 



SxMc, 



ds K^, 



• 



(31) 



It is now straightforward to read off the boundary con- 
dition by equating the coefhcient of $0 for the corre- 
sponding expression for F to zero: 

(1 • z)(^ + + i^i-w)+PK + f3Qk^0. (32) 

We have used the fact that K can be expressed as 
K = + K±, where K^^ = Kabf^t^ is the projection 
of Kab onto the unit tangent to C, In general, some 
simplification is possible by using the identity 



K\\ = n • t = cos QK\\ + sin Qk . 



(33) 



In particular, for axially symmetric geometries, we note 
that /sTii = Kn is consistent with K^^ = sin^f/i?, and 
k = cos\l//i?, with the identification 9 = + Q. 

Note that Eq. ( |3^ ) is consistent with Rcf. ||l^ where 
an axially symmetric (indeed spherical) vesicle adhering 
to a flat substrate {K = 0) is described. 



We now examine the general case, as given by the en- 
ergy (|lj) , including the bending energy Fi, . As discussed 
in Sect. Ill A, in order to avoid discontinuities at the 
boundary we impose O = 0orl-z = — lasa constraint. 

For the tangential deformation of F^, from Eq. (|2C 
we have immediately 



6i\Fb = -a / dsif^^Q. 
Jc 



(34) 



The novel non-trivial boundary term associated with 
the normal deformation of the free surface originates in 
the term 



2a f dAK-i'^H^Kab 

J free 



(35) 



contributing from the variation of F\, free. Modulo the free 
bulk shape equation (described in 122]), there remains 



h free 



2a / ds - KV 

Jc 



(36) 



We note that $ = n • z$o = 0, so that the first term 
vanishes. For the second term, on C, Eq.(p5|) gives 
V^'I' = (A'^ — A'j^)$Oj so that the novel contribution 
is 



S^_Fb fi-c 



2a / dsK [K^ - K^) $0 ■ 
Jc 



(37) 



It is now straightforward to read off the boundary con- 
dition (there is no term proportional to Vj^<f>o) 

-aK'^~(5K«+ak'^ + (5k«+2a{K^-kA_)K ^w. (38) 



ATll = when 9 = 0, 



If we now use the fact that A| 

Eq.(38) reduces to the remarkably simple expression 

\/ w/a 



Ki -A 



(39) 



This expression is independent of (3. 

The curvature K± is completely fixed at the boundary 
by Eq.(|39|). We note that the off-diagonal term with 
respect to the basis and i", A||j_ = IH^Kab = n • 1 
will not generally be zero. Just like A||, however, it will 
is completely determined by its substrate counterpart, 
A||j^ — A'li j^. Thus, all three components of the curvature 
are fixed at the boundary. 

If th e substrate is flat at the boundary, we have K±^ = 
^Jw/ a. For an axially symmetric shape, K± = Kg = 9' 
and Eq. ( p9| ) therefore reproduces the well known bound- 
ary condition |ll[. 



9' = 



(40) 



We note that if the substrate is axially symmetric and 
not flat, Eq. (EO) is modifled to 



9' - = 



(41) 



where ^' is the curvature along a meridian of the sub- 
strate. This agrees with the expression given in footnote 
14 of ITlll. 



STRESSES AT THE BOUNDARY VI. CONCLUSIONS 



In equilibrium, the forces directed along the normal 
from the boundary into the membrane must balance. In 
p2| , it is shown that the stress tensor for the model ( p^ ) 
can be expressed as 

r = [aK{2K'''' ^ K-i"^) + PiK"'' - K-i''^) 

- /i7°^]eh - 2aV°is:n. (42) 

This is the stress tensor on the free surface. It satisfies 

Va^^pn (43) 

at each point. 

Let us for simplicity suppose that the substrate is flat, 
so that Kab = 0. The corresponding stress tensor in that 
part of the vesicle which is bound to the substrate f ° is 
then 

P = («;-M)r''eb, (44) 

which is isotropic. Note that f does not satisfy the con- 
servation law Eg. (|4^) . Thus the construction of a Gaus- 
sian pillbox of infinitesimal thickness on the boundary 
does not lead to a useful identity. We note, however, 
that 1 • f°/a is the pressure acting on the boundary due 
to unbalanced stresses in the free bulk at its boundary. 
We have 

1 • f°/a = aK{2Ki_ -K)- 13K\\ - ^ . (45) 

Similarly, 

z • i^ia ^w- fi. (46) 

The stresses must balance in equilibrium. When they do, 
Eq.(^9|) is reproduced. This derivation is not only more 
efficient than the variational argument, it also homes in 
immediately on the physics at the boundary. 



We have shown how a combination of simple geomet- 
rical and variational techniques as well as conservation 
laws can be applied to study the adhesion of vesicles de- 
scribed by a geometrical hamiltonian. This provides a 
useful platform for either a numerical or perturbative ap- 
proach to adhesion, particularly when one is interested in 
non-axially symmetric shapes. Axially symmetric shapes 
are very special ones. 

These techniques also generalize to so-called floppy or 
egg-carton membranes where a term penalizing curvature 
gradients also appears in the hamiltonian ||l2|, Now, 
not only is the contact angle fixed, but its first derivative 
vanishes. It is the second derivative which will be pro- 
portional to the bond strength. The interesting shapes 
are also certainly not axially symmetric. 

Note added in proof. After this work was completed, 
Ref. was brought to our attention, where the adhe- 
sion of 'linear' vesicles in two dimensions is considered. 
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